Stellar electron-capture rates in pf-shell nuclei from quasiparticle random-phase 

approximation calculations 
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Electron-capture rates at different density and temperature conditions are evaluated for a set of pf- 
shell nuclei representative of the constituents in presupernova formations. The nuclear structure part 
of the problem is described within a quasiparticle random-phase approximation based on a deformed 
Skyrme Hartree-Fock selfconsistent mean field with pairing correlations and residual interactions in 
particle-hole and particle-particle channels. The energy distributions of the Gamow- Teller strength 
are evaluated and compared to benchmark shell-model calculations and experimental data extracted 
from charge-exchange reactions. The model dependence of the weak rates are discussed and the 
various sensitivities to both density and temperature are analyzed. 
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I. INTRODUCTION 

Stars are unique laboratories where all the interactions 
in nature come into play to determine the different stages 
in the stellar evolution. This comprises the energy gen- 
eration, which is mainly caused by fusion reactions medi- 
ated by strong interactions, as well as the nucleosynthe- 
sis of elements mediated by nuclear reactions in different 
stellar scenarios [1| ■ Whereas the main sequences of the 
stellar evolution and the majority of the elements with 
mass number below A ~ 60 are produced by hydrostatic 
burning mediated by the strong and electromagnetic in- 
teractions, weak interactions grow in importance in later 
stages |2|, when the density (p) and temperature (T) in 
the core become larger, increasing the Fermi energy of 
the degenerate electron gas and then favoring electron 
captures. 

In fact, it is a well established feature that weak /?- 
decay and electron-capture (EC) processes are very im- 
portant mechanisms to understand the late stages of the 
stellar evolution 2], playing a critical role to determine 
both the presupernova stellar structure and the nucle- 
osynthesis of heavier nuclei. These processes are domi- 
nated by Gamow- Teller (GT) transitions and therefore, 
the GT properties of pf-shell nuclei are of special im- 
portance because they are the main constituents of the 
stellar core in presupernovae formations 3] leading to 
core-collapse (type II) or thermonuclear (type la) super- 
novae. 

While the scenarios for type la supernovae are thought 
to be binary systems with a white dwarf accreting ma- 
terial from a companion star, type II supernovae are the 
final result of the gravitational collapse of the core of a 
massive star that takes place when the nuclear fuel ex- 
hausts. Then, the core becomes unstable and when the 
mass exceeds the Chandrasekhar mass, the electron de- 
generacy pressure is not able to prevent the gravitational 
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collapse. In the initial stages, electrons are captured 
by nuclei in the iron-nickel mass region, thus reducing 
the electron-to-baryon fraction (Y e ) of the presupernova 
star and correspondingly the degeneracy pressure. At the 
same time, the neutrinos produced at presupernova den- 
sities leave the star reducing the energy and cooling the 
star. Both effects act in the same direction accelerating 
the collapse. With increasing neutronization of the core 
material, the Qp energies increase and /3-decays become 
more important competing with ECs. EC processes are 
therefore essential ingredients to follow the complex dy- 
namics of core-collapse supernovae and reliable estimates 
of these rates certainly contribute to a better understand- 
ing of the explosion mechanism. 

An accurate understanding of most astrophysical pro- 
cesses requires necessarily information from nuclear 
physics, which provides the input to deal with network 
calculations and astrophysical simulations [4] . Due to the 
extreme conditions of p and T that hold in stellar scenar- 
ios, most of the nuclear properties cannot be measured 
directly. Therefore, the GT strength distributions must 
be estimated in many cases by model calculations. Obvi- 
ously, nuclear physics uncertainties will finally affect the 
reliability of the description of those astrophysical pro- 
cesses. 

The first extensive calculations of stellar weak rates as 
functions of relevant ranges of p and T were done in Ref. 
[2j . It was assumed that the whole GT strength resides in 
a single resonance whose energy relative to the daughter 
ground state is parametrized phenomenologically, taking 
the total GT strength from the single-particle model. 

In the last decades, GT + strength distributions on nu- 
clei in the mass region A ~ 60 have been studied ex- 
perimentally via (n,p), or equivalent higher resolution 
charge-exchange reactions such as (d, 2 He) and (i, 3 He), 
at forward angles [5H20J. Charge- exchange reactions are 
the most efficient way to extract the GT + strength in sta- 
ble nuclei |21|. For incident energies above 100 MeV, the 
isovector spin-flip component of the effective interaction 
is dominant and the cross sections is mainly originated 
from spin-isospin transitions. At forward angles the mo- 



mentum transfer is small and the reaction cross section 
is dominated by the GT operator with AT = 1,AL = 
0, A J w = 1 + . The cross section, extrapolated to zero mo- 
mentum transfer, is proportional to the /3-decay strength 
between the same states. Charge-exchange reactions at 
small momentum transfer are therefore used to study GT 
strength distributions when /3-decay is not energetically 
possible. 

The data show that the total GT + strength is strongly 
quenched and fragmented over many final states, as com- 
pared to the independent-particle model. Then, improve- 
ments in the weak rates have been focused on the de- 
scription of the nuclear structure aspect of the problem. 
Different approaches to describe the nuclear structure 
involved in the stellar weak decay rates can be found 
in the literature. They are basically divided into shell 
model (SM), either SM Montecarlo [22[ or large scale SM 
diagonalizations [23-25] and proton-neutron quasiparti- 
cle random-phase approximation (QRPA) [26-31] cate- 
gories. Although QRPA calculations cannot reach the 
detailed spectroscopy achieved from present state-of-the- 
art SM calculations, the global performance of QRPA 
is quite satisfactory. Moreover, one clear advantage of 
the QRPA method is that it can be extended to heav- 
ier nuclei, which are beyond the present capability of full 
SM calculations, without increasing the complexity of the 
calculation. 

Very recently, a systematic evaluation of the ability 
to reproduce the measured GT strength distributions of 
various theoretical models based on SM and QRPA was 
done in Ref. [32j, where EC rates were derived from 
those models at relevant p and T. While several sets 
of SM calculations using different effective interactions 
were compared, namely KB3G [33| and GXPFla [3j, in 
the case of QRPA only the formalism developed in 



using deformations and masses obtained from the finite 
range droplet model [36[ was considered in Ref. [33|. In 
what follows we use the term QRPA-Mdller to refer to 
these QRPA results obtained from the above mentioned 
formalism. 

Given the sensitivity of the weak rates to the nuclear 
structure through the GT strength distributions, it is 
worth extending the study in Ref. [32] by considering 
alternative calculations within the QRPA formalism. In- 
deed, the QRPA method with separable GT interactions 
was first proposed and applied in Ref. [371 ] . on a spher- 
ical harmonic oscillator basis, and then it was extended 
to deformed nuclei |35j using deformed phenomenological 
single-particle basis. Further refinements to the QRPA 
formalism were introduced along the years [38l - l44fl . 

In this work, we study the dependence of the EC rates 
on both p and T with GT strength distributions calcu- 
lated within a QRPA approach based on a selfconsis- 
tent deformed Hartree-Fock (HF) mean field with Skyrme 
interactions including pairing correlations and residual 
separable forces in both particle-hole (ph) and particle- 
particle (pp) channels. We compare our calculations with 
the benchmark calculations in Ref. [331. This formalism 



represents an improvement over the QRPA-Moller ap- 
proach in several aspects. First, instead of a phenomeno- 
logical approach (single-particle models based on Nilsson, 
Woods-Saxon, or folded Yukawa models), the deformed 
mean field is now obtained selfconsistently and there is no 
need to introduce deformation parameters as input, and 
secondly, a separable residual GT interaction in the pp 
channel is included. We should also mention the inclusion 
of an effective quenching factor in our calculations that 
was not included in QRPA-Moller. The present nuclear 
model has been tested successfully reproducing very rea- 
sonably the experimental information available on both 
bulk and decay properties of medium- mass nuclei [45l450| . 
The paper is organized as follows. In Section [TT] the 
weak decay rates are introduced as functions of p and T 
and their nuclear structure and phase space components 
are described. Section IIIII contains the results obtained 
for the GT strength distributions and for the EC rates 
in some selected p/-shell nuclei that correspond to those 
used in Ref. [32[ . Section IIVI contains the conclusions of 
this work. 



II. THEORETICAL FORMALISM 

There are several distinctions between terrestrial and 
stellar decay rates caused by the effect of high p and T. 
One effect of T is directly related to the thermal popu- 
lation of excited states in the decaying nucleus, accom- 
panied by the corresponding depopulation of the ground 
states. The weak-decay rates of excited states can be 
significantly different from those of the ground state and 
a case by case consideration is needed. Another effect 
related to the high p and T comes from the fact that 
atoms in these scenarios are completely ionized and con- 
sequently electrons are no longer bound to the nuclei, but 
forming a degenerate plasma obeying a Fermi-Dirac dis- 
tribution. This opens the possibility for continuum EC, 
in contrast to the orbital EC caused by bound electrons 
in the atom under terrestrial conditions. These effects 
make weak interaction rates in the stellar interior sensi- 
tive functions of T and p. 



A. Weak decay rates 

Assuming thermal equilibrium, the probability of occu- 
pation of the excited states in the parent nucleus follows 
a Boltzmann distribution. The decay rate of the parent 
nucleus is given by 



A = 5> 
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where G = Y\ (2«7j + 1) e - B °/( fc s T ) is the partition func- 
tion and Ji(Ei) is the angular momentum (excitation en- 
ergy) of the parent nucleus state i. 



In principle, the sum extends over all populated states 
in the parent nucleus up to the proton separation energy. 
However, because of the range of temperatures consid- 
ered in this work (T = 1 — 10 GK), only a few low- lying 
excited states are expected to contribute in the decay 
of even-even nuclei. In fact, the lowest excited states in 
the even-even nuclei considered are collective 2 + states 
located typically above 1 MeV from the ground state. 
The lowest of these 2 + states corresponds to the nucleus 
56 Fe and appears at E 2 + = 0.847 MeV [5l|. Hence, their 
contributions to the rates can be neglected at these tem- 
peratures. 

In the case of odd-A nuclei the situation is more in- 
volved because excited states of single-particle nature ap- 
pear at low excitation energies. This is particularly the 
case of 45 Sc, where one finds up to six excited states be- 
low 1 MeV [51j. These states will become populated as 
T raises and may contribute to the weak rates. A case 
by case analysis is mandatory and work in this line is in 
progress. 

The decay rate corresponding to the parent state i is 
given by 
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where the sum extends over all the states in the final 
nucleus reached in the decay process and D = 6146 s. 
This expression is decomposed into a phase space factor 
<&if, which is a function of p and T and a nuclear struc- 
ture part Bif that contains the transition probabilities 
for allowed Fermi and GT transitions, 



B if = B if (GT)+B if (F). 



(3) 



In this work we only consider the dominant GT tran- 
sitions. Fermi transitions have a simple expression as- 
suming isospin symmetry and are only important for j3 + 
decay of neutron-deficient light nuclei with Z > N. The 
theoretical description of both Bif and $i/ are explained 
in the next subsections. 



with Tg = 10 9 K) Contributions to the weak rates from 
positron decays are therefore neglected in this work. 

In the astrophysical scenarios of our study, nuclei are 
fully ionized and continuum EC from the degenerate elec- 
tron plasma are possible. The phase space factor for EC 
is given by 



$f/ = / MQ^+LofFiZ^) 

xS e {uo)[l-S v {Q lf +uj)]duj. (4) 

In this expression t o is the total energy of the electron in 
m e c 2 units, p = yup — 1 is the momentum, and Qif is 
the total energy available in m e c 2 units 
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written in terms of the nuclear masses of parent (M p ) 
and daughter (M^) nuclei and their excitation energies 
E t and Et, respectively. 

F(Z, uj) is the Fermi function that takes into account 
the distortion of the electron wave function due to the 
Coulomb interaction. 



F(Z,u) = 2(1 + 7 )(2pi?)- 2 ( 1 ^)e^fe±MjJ , (7) 



where 7 = yl — (qZ)' 
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aZuj /p ; a is the fine 



structure constant and R the nuclear radius. The lower 
integration limit in Eq. Q is given by cue = 1 if Qif > 
-1, or uj ( = \Q if \ if Qif < -1. 

S e and S v , are the electron and neutrino distribution 
functions, respectively. Its presence inhibits or enhances 
the phase space available. In the stellar scenarios con- 
sidered here the commonly accepted assumption is that 
S v = 0, because neutrinos and antineutrinos at these 
densities can escape freely from the interior of the star. 
The electron distribution is described as a Fermi-Dirac 
distribution 



B. Phase Space Factors 

The nuclei under study in this work correspond to sta- 
ble p/-shell nuclei. f3 + decays from their ground states 
are then energetically forbidden. However, as T raises, 
thermal population of excited states in the parent nu- 
cleus may induce /3 + -decays if Ei exceeds the Qa+ en- 
ergy. These decays that are almost independent of p 
and T, might compete with ECs in cases where T is high 
enough to populate significantly excited states in the par- 
ent nucleus. In our case the decays would involve excited 
states beyond 1.5-2 MeV in the most favored cases, but 
they are not sufficiently populated even at the higher T 
considered in this work (k B T = 0.862 MeV at T 9 = 10 
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exp [(w - Me) /{k B T)\ + 1 



(8) 



The chemical potential fi e as a function of p and T is 
determined from the expression 



pye= ^k(^) 3 I p--w*> 



(9) 



in (mol/cm 3 ) units, p is the baryon density (g/cm 3 ), Y e 
is the electron-to-baryon ratio (mol/g), and Na is Avo- 
gadro's number (mol -1 ). 

The phase space factor for EC in Eq. (UJ) is therefore a 
sensitive function of both p and T, through the electron 
distribution S e . 



C. Nuclear Structure 

The nuclear structure part of the problem is described 
within the QRPA formalism. Various approaches have 
been developed in the past to describe the spin-isospin 
nuclear excitations in QRPA [29 l l35l. I38U44J . In this sub- 
section we show briefly the theoretical framework used 
in this work to describe the nuclear part of the decay 
rates. More details of the formalism can be found in 
Refs. (SSI S3. 

The method starts with a self-consistent deformed 
Hartree-Fock mean field formalism obtained with Skyrme 
interactions, including pairing correlations. The single- 
particle energies, wave functions, and occupation proba- 
bilities are generated from this mean field. In this work 
we have chosen the Skyrme force SLy4 [52| as a represen- 
tative of the Skyrme forces. It is one of the most success- 
ful Skyrme forces and has been extensively studied in the 
last years. We also consider for comparison the results 
obtained with the force SG2 (53| that has been success- 
fully tested against spin-isospin excitations in spherical 
and deformed nuclei. 

The solution of the HF equation is found by using the 
formalism developed in Ref. [54J , assuming time reversal 
and axial symmetry. The single-particle wave functions 
are expanded in terms of the eigenstates of an axially 
symmetric harmonic oscillator in cylindrical coordinates, 
using twelve major shells. The method also includes pair- 
ing between like nucleons in BCS approximation with 
fixed gap parameters for protons and neutrons, which are 
determined phcnomenologically from the odd-even mass 
differences involving the experimental binding energies 
51]. 

Potential energy curves are analyzed as a function of 
the quadrupole deformation. For that purpose, con- 
strained HF calculations are performed with a quadratic 
constraint [55]. The HF energy is minimized under the 
constraint of keeping fixed the nuclear deformation. Cal- 
culations for GT strengths are performed subsequently 
for the various equilibrium shapes of each nucleus, that 
is, for the solutions, in general deformed, for which min- 
ima are obtained in the energy curves. 

To describe GT transitions, a spin-isospin residual in- 
teraction is added to the Skyrme mean field and treated 
in a deformed proton-neutron QRPA. Taking into ac- 
count the phenomenological nature of the Skyrme inter- 
actions, which are not fitted to spin-isospin observables, 
it is reasonable to complement the Skyrme interaction 
with some extra phenomenological parameters sensitive 
to those nuclear properties. 

This interaction contains a ph and a pp part. The in- 
teraction in the ph channel is mainly responsible for the 
position and structure of the GT resonance and it can 
be derived consistently from the same Skyrme interac- 
tion used to generate the mean field, through the second 
derivatives of the energy density functional with respect 
to the one-body densities. The ph residual interaction 
is finally expressed in a separable form by averaging the 



resulting contact force over the nuclear volume [40(. By 
taking separable GT forces, the energy eigenvalue prob- 
lem reduces to find the roots of an algebraic equation. 

The pp part is a neutron-proton pairing force in the 
J^ = 1 + coupling channel, which is also introduced as a 
separable force [38|, |42| . Its strength is usually fitted to 
reproduce globally the experimental half-lives. Various 
attempts have been done in the past to fix this strength 
[5a], arriving to expressions that depend on the model 
used to describe the mean field, Nilsson model in the 
above reference. In previous works J40h431 |45j we have 
studied the sensitivity of the GT strength distributions 
to the various ingredients contributing to the deformed 
QRPA calculations, namely to the nucleon-nucleon effec- 
tive force, to pairing correlations, and to residual inter- 
actions. We found different sensitivities to them. In this 
work, all of these ingredients have been fixed to the most 
reasonable choices found previously. In particular we use 
the coupling strengths Xgt = 0.10 MeV and k P q T = 0.05 
MeV. An optimum set of coupling strengths (xgt> kP gt) 
could be chosen following a case by case fitting proce- 
dure and we will get different answers depending on the 
nucleus, shape, and Skyrme force. However, since the 
purpose here is to test the ability of QRPA to account 
for the GT strength distributions in the iron-nickel mass 
region with as few free parameters as possible, we have 
chosen to use the same coupling strengths for all the nu- 
clei considered in this work. 

The GT strength for a transition from an initial state 
i to a final state / is given by 



Bij(GT±) = 



2J t 
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(10) 



(gA/gv)ctt = Q -7 {g a / gv)ba.re is the effective ratio of axial 
and vector coupling factors that takes into account in 
an effective manner the observed quenching of the GT 
strength. 

In even-even nuclei, the GT strength corresponding to 
the transition JiKi — > JfKf (0 + — > \ + K) in the labora- 
tory system, is expressed in terms of the intrinsic ampli- 
tudes connecting the QRPA ground state \4>$) with the 
one-phonon excited states \4>k) 



B if (GT±) 



9v/ cS 



8k,o(<Pk Iffo^l^o) (11) 
+2 6 KA {^ K \a 1 t ± \ ( j )Q ) 2 }{12) 



To obtain this expression, the initial and final states in 
the laboratory frame have been expressed in terms of 
the intrinsic states using the Bohr-Mottelson factoriza- 
tion [13] ■ 

When the parent nucleus has an odd nucleon, the 
ground state can be expressed as a one quasi-particle 
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FIG. 1: (Color online) Gamow- Teller strength distribution 
B(GT+) for the transition 45 Sc to 45 Ca plotted versus the 
excitation energy of the daughter nucleus, (a) Experimental 
data are compared to our QRPA results with SLy4 Skyrme 
force (individual transitions) as well as with folded distribu- 
tions from SLy4 and SG2 forces, (b) Experimental accu- 
mulated B(GT + ) strength compared to QRPA results from 
Moller [35|] as given in Ref. [32] and from our calculations 
with SLy4 and SG2 interactions, (c) Experimental accumu- 
lated B(GT + ) strength compared to shell-model results with 
KB3G and GXPFla interactions. Data extracted from (n,p) 
reactions are from Ref. Q. 



state in which the odd nucleon occupies the single- 
particle orbit of lowest energy. Then two types of tran- 
sitions are possible. One type is due to phonon excita- 
tions in which the odd nucleon acts only as a spectator. 
In the intrinsic frame, the transition amplitudes are in 
this case basically the same as in the even-even case but 
with the blocked spectator excluded from the calcula- 
tion. The other type of transitions are those involving 
the odd nucleon state, which are treated by taking into 
account phonon correlations in the quasiparticle transi- 
tions in first order perturbation J35J, |38|, |42| . 



III. RESULTS 

A. Gamow- Teller distributions 

In the next figures, we show the results obtained for the 
energy distributions of the GT strength corresponding to 
the equilibrium shapes for which we obtained minima in 



FIG. 2: (Color online) Same as in Fig. [TJ but for the 48 Ti 
to 48 Sc transition. Data are from (n,p) [6j] and (d, 2 He) Q 
reactions. 
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FIG. 3: (Color online) Same as in Fig. |T] but for the 51 V 
to jl Ti transition. Data are from (n,p) |3] and (d, He) [§] 
reactions. 
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FIG. 4: (Color online) Same as in Fig. [T] but for the 54 Fe 
to 54 Mn transition. Data are from (n,p) reactions from Refs. 
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FIG. 6: (Color online) Same as in Fig. [T] but for the 56 Fe 
to 56 Mn transition. Data are from (n,p) reactions from Refs. 
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FIG. 5: (Color online) Same as in Fig. Q] but for the 55 Mn 
to 55 Cr transition. Data are from (n,p) reactions [12 ]. 
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FIG. 7: (Color online) Same as in Fig. Q] but for the 58 Ni to 
58 Co transition. Data are from (n,p) [12], (d, 2 He) [1^]. and 
(t, 3 He) [3] reactions. 
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FIG. 8: (Color online) Same as in Fig. Q] but for the 59 Co to 
59 Fe transition. Data are from (n,p) reactions |8[. 



FIG. 10: (Color online) Same as in Fig. [T] but for the 62 Ni 
to 62 Co transition. Data are from (n,p) reactions [15| . 
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FIG. 9: (Color online) Same as in Fig. [Q but for the 60 Ni 



to 
15 



Co transition. Data are from (n,p) reactions from Refs. FI ^ H= (Color online) Same as in Fig. m but for the 64 m 
rai to Co transition. Data are from (n, p) [15(] and (d, He) [171 ] 

reactions. 
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FIG. 12: (Color online) Same as in Fig. [TJ but for the 64 Zn 
to 64 Cu transition. Data are from (d, 2 He) [lg| and (£, 3 He) 
[igj ] reactions. 



the potential energy curves. The GT strength is plot- 
ted versus the excitation energy of the daughter nucleus 
E ex = E f (MeV). 

Figs. [TlfTSI contain the results for the isotopes 45 Sc, 
48 Ti, 51 V, 54 Fe, 55 Mn, 56 Fe, 58 Ni, 59 Co, 60 Ni, 62 Ni, 64 Ni, 
and 64 Zn, respectively. This set of nuclei, which are the 
object of our study, have been chosen in accordance with 
the cases studied in Ref. |32j |. 

In the top panels (a) we show the experimental en- 
ergy distributions and the individual GT strengths cor- 
responding to the SLy4 interaction. We also show the 
continuous distributions from SLy4 and SG2 obtained 
by folding the strength with Breit-Wigner functions as it 
was done in Refs. [4^, [58, 59], so that the original discrete 
spectrum is transformed into a continuous profile. In the 
middle panels (b) we plot the measured energy distribu- 
tion of the accumulated strength, which is compared to 
three QRPA calculations. They are the QRPA-Moller 
calculations and the two calculations obtained from our 
present formalism using the Skyrme forces SLy4 and 
SG2. Finally, in the bottom panels (c) the comparison of 
the experimental distribution is made with SM calcula- 
tions with the interactions KB3G [H and GXPFla fil . 
as given in Ref. 



32|. 



Fig. Q] shows the GT strength distributions for the 
45 Sc(7/2~) — > 45 Ca extracted from (n,p) charge-exchange 
reactions [5J. Results are available in 1 MeV energy bins. 
The measured strength distribution is compared to our 



QRPA results in Fig. [lja). The structure of the calcu- 
lated distribution agrees qualitatively with the data in 
the sense that the GT strength is concentrated around 
7 MeV with small bumps at low energy. Both SG2 and 
SLy4 results are very similar. The total accumulated 
strength as a function of the excitation energy can be 
seen in Fig. |TJb) , where the running sums are compared. 
Our calculations follow closely the experimental strength, 
while QRPA-Moller overestimates it. Similarly, the SM 
calculations shown in Fig. QJc) compare reasonably well 
with the experiment with the strength practically concen- 
trated in a single peak at 6 MeV and somewhat smaller 
total strength in the whole energy range. 

Fig. [2] contains the results for the 48 Ti(0+) -> 48 Sc 
transition. In this case the data are from (n,p) [6j and 
(d 2 He) [3] reactions. One should note that the strength 
measured in Ref. [6j would contain an isovector spin 
monopole component, which is estimated to contribute 
to about one third of the total strength. In the case of 
the (d 2 He) reaction, the high resolution achieved was 
120 keV, but the strength was extracted only up to an 
excitation energy of 5 MeV. Both sets of data agree in 
the location of the peak of the strength distribution at 
about 3 MeV, but disagree in other regions. Our QRPA 
calculations show a two bump structure with peaks cen- 
tered at about 4 and 6 MeV that resembles the profile of 
the (n,p) data. The total strength obtained lies between 
the strengths from (n,p) and (d 2 He) reactions and is 
lower than QRPA-Moller. The SM calculations accumu- 
late the strength around 4 MeV containing less strength 
than QRPA. 

Fig. E] for the 51 V(7/2~) -^ 51 Ti transition contains 
data from (n,p) [8| and (d 2 He) reactions with 1 MeV 
and 110 keV resolution, respectively. Strengths from the 
high resolution experiments were only extracted up to 6.5 
MeV. The main characteristic of these data is that the 
strength appears concentrated at about 5 MeV excitation 
energies, a feature that is well reproduced by our QRPA 
results although at somewhat higher energy. The total 
strength measured in both experiments are similar and 
agrees quite well with our QRPA and with SM results. 

The results for the 54 Fe(0 + ) -> 54 Mn transition appear 
in Fig. |U In this case two sets of data from (n,p) re- 
actions are available from Refs. [RJ, LU|- Our QRPA 
results produce a bump centered at about 5 MeV which 
is a little bit displaced to higher energies with respect to 
the experiment. The total strengths in QRPA are also 
above the measured strength, which is better reproduced 
by the SM calculations. 

In the case of the 55 Mn(5/2~) — > 55 Cr transition in Fig. 
[5] there are data extracted from (n,p) reactions [l2|. The 
broad peak observed experimentally centered at about 4 
MeV is well reproduced in our QRPA calculations. Also 
the total GT strength is in this case well accounted for 
by all the models. 

Fig. [5] shows the results corresponding to the 
56 Fe(0 + ) — » 56 Mn transition. Data are from (n,p) re- 
actions do, nJI . Both sets of data show a concentration 



of the GT strength between 1 and 4 MeV in agreement 
with the calculations. However, the strength observed 
beyond 6 MeV is not found in any of the calculations 
presented. The total strength measured is slightly over- 
estimated (underestimated) by QRPA (SM) calculations. 

The distribution of the GT strength corresponding to 
the transition 58 _Ni(0 + ) -^ 58 Co is shown in Fig. [71 Data 
are from (n,p) [1 2:] , (d 2 He) [llij . and (£, 3 He) [l4( reac- 
tions with energy resolutions of 1.2 MeV, 130 keV, and 
250 keV, respectively. Our QRPA calculations produce 
a strength distribution sharply concentrated between 3 
and 4 MeV that contains practically all the GT strength 
observed. The SM results are more fragmented in better 
agreement with experiment. Contrary to previous cases, 
we use here and in the rest of Ni isotopes, the SM results 
obtained with the improved interaction GXPF1J [25J. 

Fig. |5] shows the results for the 59 Co(7/2~) -> 59 Fe 
transition. Data in this case have been obtained from 
(n,p) reactions [8] with about 1 MeV energy resolution. 
QRPA results show one peak structure centered around 
5 MeV, which seems to be displaced 1 MeV to higher en- 
ergy with respect to the experimental distribution. The 
total GT strength is somewhat larger than experiment. 
On the other hand SM results agree quite well with the 
experiment. 

In the case of the 60 Ni(0 + ) — ^ 60 Co transition shown 
in Fig. [SI the data are from (n,p) reactions [15| and 
from the reanalysis performed in Ref. [la ]. As in the 
case of 58 Ni, whereas the experimental strength appears 
fragmented below 3 MeV, QRPA calculations produce 
a sharp transition with practically all the strength at 2 
MeV. All the calculations except SM-KB3G overestimate 
the total strength measured. 

The case of the 62 Ni(0+) -^ 62 Co transition (Fig. [TDJ is 
very similar to the previous case for 60 Ni. The data from 
(n,p) reactions 15] show that practically all the strength 
is contained below 2 MeV, which is compatible with the 
structure of the GT strength distributions obtained from 
the theoretical models. However, while QRPA calcula- 
tions overestimate the total strength, SM results agree 
very nicely with experiment. 

For the 64 Ni(0+) -^ 64 Co transition in Fig. [TT]we have 
data not only from (n,p) reactions [15J as in the previous 
Ni isotopes, but also data up to 4 MeV from (d 2 He) (l7| 
reactions at a much higher energy resolution (110 keV). 
Most of the GT strength is observed in the ground state 
to ground state transition, while in QRPA the strength 
is fragmented below 2 MeV. The calculated strength is 
overestimated (underestimated) by QRPA (SM) calcula- 
tions. 

In the last example, we see in Fig. [I2]the results for the 
64 Zn(0+) ^ 64 Cu transition. Data are from (d 2 He) pj] 
and (i, 3 He) [19| reactions achieving 115 keV and 280 keV 
energy resolution, respectively. The QRPA results show 
a bump centered at an excitation energy of 3 MeV (5 
MeV) when using SG2 (SLy4) interaction. SG2 results 
reproduce better the experimental results, which show 
a strong peak at 3 MeV. The calculated total strength, 



including results from SM calculations, overestimates al- 
ways the experiment. Results from SM-KB3G produce a 
strong peak at 5 MeV at variance with experiment. 

Summarizing this section we can say that the observed 
fragmentation of the experimental GT + strength distri- 
butions over many states, as well as the centroids and 
widths of the distributions are reasonably well described 
by the calculations. In the experimental distributions, 
especially in those extracted from (n,p) reactions, there 
is a tendency to build up a second peak beyond ~ 6 MeV, 
which is not reproduced in the calculations. In general, 
QRPA produce more strength at higher energy than SM 
because of the higher N-shell mixing contained in QRPA. 
All in all, the present QRPA calculations based on the 
deformed Skyrme HF+BCS+QRPA described in Sec. [H] 
are in general of comparable quality to SM calculations. 
These results serve to refine the findings in Ref. p32j |. 
where it was concluded that QRPA calculations based 
on Ref. [35j produce systematically much larger devia- 
tions from the data than SM calculations. We have shown 
in this work that this is not necessarily the case for all 
QRPA type calculations. 
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FIG. 13: (Color online) Electron-capture rates for 45 Sc 
obtained from the experimental GT strength distributions 
and from different shell model and QRPA calculations as 
a function of the temperature Tg (GK) for densities pY c — 
10 10 , IO 9 , IO 8 , IO 7 , IO 6 mol/cm 3 in panels (a), (b), (c), (d), 
and (e), respectively. 



Stellar weak decay rates 



FIG. 14: (Color online) Same as in Fig. [TJ1 but for 48 Ti. 



In the following figures (Figs. [I~3ll24]) we present the 
EC rates of the selected p/-shell nuclei as a function of 
the temperature and for various densities. The range of T 
considered varies from Tg = 1 up to Tg — 10, whereas the 
range in pY e varies from pY e = IO 6 mol/cm 3 up to pY e = 
10 10 mol/cm 3 . This grid of p and T includes those ranges 
relevant for astrophysical scenarios related to the silicon- 
burning stage in a presupernovae star [3[ (pY e = IO 7 
mol/cm 3 and Tg = 3), as well as scenarios related to pre- 
collapse of the core [60J and thermonuclear runaway type 
la supernovae [6l| {pY e = 10 9 mol/cm 3 and Tg = 10). 
In each figure from Fig. [13] to Fig. [24] we show the 
EC rates obtained from the experimental GT strength 
distributions and from different SM (KB3G and GXPF1) 
and QRPA (SLy4, SG2, and Moller) calculations. 

Concerning the experimental EC rates AEC,exp in Figs. 
n"3l!Hl one should keep in mind that these quantities are 
not necessarily the actual rates in stellar scenarios at high 
p and T. Aec.cxp is indeed the rate calculated from the 
experimental GT strength distribution extracted from 
charge-exchange reactions. Besides the intrinsic uncer- 
tainties in the extraction of the GT strength due to sev- 
eral causes like the global normalization of the unit cross 
section or to possible interferences caused by the tensor 
component of the interaction, which are typically esti- 
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FIG. 15: (Color online) Same as in Fig. \M but for 51 V. 



mated to be a 10-20% effect, one has to consider other 
sources of ignorance that makes Aec.cxp to be different 
from the EC rates in stellar scenarios. First, the GT 
strength is only measured up to some excitation energy 
and therefore, AEC.exp does not include possible contribu- 
tions from transitions beyond the measured energy range 
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FIG. 16: (Color online) Same as in Fig. \M but for 54 Fe. 



that could have an effect, especially at high T and p. 
Secondly, even when the GT strength distribution is per- 
fectly determined from charge-exchange reactions in the 
laboratory under terrestrial conditions, this is not suffi- 
cient to determine the EC rates in stellar scenarios where 
ECs can occur on excited states of the parent nucleus 
that become thermally populated at high T. Although 
these effects are expected to be very small for the den- 
sities and temperatures considered, one should keep in 
mind the real significance of Aec,cxp- 

Some general comments about the sensitivity of the EC 
rates to (p,T), to the Fermi and Qec energies, and to 
the GT distribution are in order to understand better the 
EC rates. Values for Qec and Fermi energies are given 
in Table I and Table II, respectively. Since the Fermi 
energy increases with the density, it is expected that the 
EC rates at low densities are mainly sensitive to the GT 
strength of states at low excitation energies in the daugh- 
ter nucleus. This is especially true at low T, where the 
shape of the electron energy distribution, S e , has a sharp 
surface at the Fermi energy. When T raises, the shape 
of S e is smeared out and thus, GT transitions at higher 
excitation energies might contribute to the rates even at 
low densities. When p increases the rates also increase 
because the Fermi energy of the electrons is larger and 
larger allowing to reach higher excitation energies in the 
daughter nucleus and thus making the GT strength at 
these energies contribute to the rates in a more signifi- 
cant way. The magnitude of this general effect is different 
in each nucleus because of the different Qec values that 
force the electron energy to be large enough to overcome 
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FIG. 17: (Color online) Same as in Fig. [TJ but for 55 Mn. 
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FIG. 18: (Color online) Same as in Fig. \M but for 56 Fe. 



it. When increasing T, the rates in general increase be- 
cause the diffuseness of the electron energy distributions 
makes it possible that more excited states at higher en- 
ergy can be reached. 

On the other hand the role of the Qif energy in Eq. ([5]) 
is also very important. It determines the lower integra- 
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FIG. 19: (Color online) Same as in Fig. \M but for 58 Ni. 



tion limit in the EC phase factor |QJ. Then, the energy 
of the available electrons must overcome this value to be 
captured or in other words, the Fermi energy has to be 
larger than Qij. The argument is strictly valid at T = 
where S e changes abruptly from one to zero at the Fermi 
energy. When T increases, S e is smeared out and ECs 
are possible even for Fermi energies lower than Qij. In 
general, the larger the Qec energy (more negative) the 
lower the EC rates. This effect will be accentuated at 
low densities where the Fermi energy is small. Similarly, 
nuclei with large (negative) Qec values (i.e., 48 Ti, 56 Fe, 
62 Ni, 64 Ni) are mostly sensitive to the GT strength of the 
ground and lowest excited states and the opposite is true 
for the nuclei with small Qec values ( 45 Sc, 54 Fe, 58 Ni, 
64 Zn). 

It is also worth comparing the rates in the Ni isotopes 
that exhibit different values of their Qec energies (see 
Table I). The EC rates in 64 Ni {Q EC = -7.307 MeV) 
are almost independent of the density for values between 
10 6 and 10 8 . They start growing at 10 9 and become very 
large in comparison at 10 10 . This effect is related to the 
fact that the Fermi energy is about 11 MeV (see Table 
II) at the highest density, which is enough to surpass the 
large Qec value. We find a similar situation in the case 
of 62 Ni {Qec = -5.315 MeV) and to a less extent in 60 Ni 
(Q EC = -2.823 MeV). Finally, in 58 Ni (Q EC = -0.382 
MeV) the rates increase more steadily because the Fermi 
energy exceeds the small Qec much sooner. Nuclei with 
large negative Qec energies such as 48 Ti and 56 Fe show 
similarities with the rates in 64 Ni, whereas nuclei with 
small negative Qec energies like 45 Sc, 54 Fe, and 64 Zn 



FIG. 20: (Color online) Same as in Fig. [H but for 59 Co. 
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FIG. 21: (Color online) Same as in Fig. OH but for 60 Ni. 



exhibit similarities with the rates in 58 Ni. 

We should also mention that in all the examples stud- 
ied, the rates are practically independent of T at pY e = 
10 10 mol/cm 3 and only a model dependence is apparent. 
This is a direct consequence of the large Fermi energy 
(about 11 MeV) at any T that makes the rates sensitive 
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FIG. 22: (Color online) Same as in Fig. [H but for 62 Ni. 
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FIG. 23: (Color online) Same as in Fig. H3 but for 64 Ni. 



to the GT strength at all the excitation energies. 

Now we discuss the model dependence of the EC rates. 
As we have already mentioned, as a general rule the EC 
rates at low p and T will be more sensitive to the GT 
strength at low excitation energies, especially for nuclei 
with large negative Qec energies. On the other hand, 




FIG. 24: (Color online) Same as in Fig. [TJ but for 64 Zn. 



the rates at high p and T will depend more on the global 
structure of the GT distribution. This type of correla- 
tions can be seen to some extent in all the figures from 
Figs [TJ to Fig. [Ml 

In the case of 45 Sc (Fig. Il~3l) the small value of Qec 
makes the EC rates very sensitive to the nuclear model 
because the GT distribution is to a large extent involved 
up to large excitation energies. We can see that we always 
get the rates from QRPA-Moller larger than the rates 
from QRPA-SLy4 and QRPA-SG2, and larger than the 
SM ones. This is clearly correlated to the GT strength 
distributions in Fig. [T] where we observe the same or- 
dering in the GT strength. In particular, the rates from 
our QRPA calculations are very close to the experimen- 
tal ones, while QRPA-Moller gives larger rates and SM 
produce rates clearly above. This is especially noticeable 
at low p and T due to the absence of strength between 
0-6 MeV in the SM calculations. 

In the case of 48 Ti (Fig. Q3J) with a relatively large 
Qec value, we can see that the experimental rates from 
(n,p) are larger than the corresponding rates for (d, 2 He) 
in accordance with the GT strengths in Fig. [2] While 
the QRPA calculations agree better with the rates from 
(n,p), the SM calculations agree better with the rates 
from (d, 2 He). In the case of 51 V (Fig. ITS]) all the models 
are quite similar with the exception of QRPA-Moller that 
produces larger rates, a feature that is connected with the 
structure of the GT strength distribution in Fig. ([3]). 

For 54 Fe (Fig. [TB")) Qec is again very small allowing 
most of the GT strength to be probed more easily. At low 
densities, the agreement of our QRPA calculations with 



TABLE I: Qec values [MeV] from Eq. (J6| using experimental nuclear masses from Ref. [51 

45g c 48 Ti Sly 54p e 55 Mn 56p e 58 M 59^ 60 Ni BBjjj 6* Ni 64^ 

-0.258 -3.990 -2.472 -0.697 -2.603 -3.696 -0.382 -1.565 -2.823 -5.315 -7.307 -0.580 
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experiment is better than for other models. It is also 
interesting to observe how the rates from QRPA-Moller 
are much lower than other models at low p and T, but 
are larger at the largest density. This is a consequence of 
the GT distribution in Fig. [4] that shows a huge strength 
at 4 MeV and practically nothing else. Then, at low p 
and T there are no electrons available to be captured to 
that energy. On the contrary at high densities, when the 
Fermi energy is larger, strong EC to that state is possible 
increasing dramatically the rate. 

In the case of 55 Mn (Fig. ITTf the most remarkable 
aspect to mention is the similarity between the model 
predictions, as well as the discrepancy with experiment 
at low p and T. This is caused by the larger experi- 
mental GT strength at low energies. At higher p when 
the effect of the low-lying states is not so important the 
agreement with the experiment improves. The situation 
is very similar in 56 Fe (Fig. IT51) . where the experimental 
GT strength at low excitation energy is larger than any 
model prediction. The rates from QRPA-Moller are par- 
ticularly low as a consequence of the lack of GT strength 
below 2.5 MeV (see Fig. |SJ}. 

In 58 Ni (Fig. [H))) with a very small Qec we can see 
that at low p the rates are sensitive to the strength of 
the low-lying excitations and thus the rates from QRPA- 
Moller are the smallest, while at the high p the rates are 
sensitive to the total strength and then QRPA-Moller is 
the largest in accordance with the GT strengths in Fig. [7] 
The situation in 59 Co (Fig. [20|) is similar to the previous 
case, but in this case our QRPA rates are the lowest at 
low p because the GT strengths are the smallest at low 
energies. 

In the cases of 64 Ni and 62 Ni, and to a less extent 60 Ni, 
Qec has large negative values and therefore the rates are 
mainly sensitive to the strength of the low-lying states. 
In the case of 60 Ni (Fig. |2"TJ) the rates from the different 
models at low p vary from the larger values of the exper- 
imental rates up to the lower values of the QRPA-Moller 
calculations. This is in agreement with the GT strength 
at low energies, where the measured strength is larger 
than any model. The GT strength from QRPA-Moller is 
practically inexistent up to 2 MeV. At higher p one finds 
the rates ordered from the higher QRPA to the lower SM, 
passing through the experimental rates. In this case this 
ordering follows the total GT strength in Fig. [9] as one 
could expect from the fact that at high p, the Fermi en- 
ergy is high enough to overcome the range of excitation 
energies where GT strength is found. Similar arguments 




FIG. 25: (Color online) Ratios of various EC rates calculated 
from theoretical and experimental GT strength distributions 
with respect to the EC rates calculated in QRPA with the 
Skyrme force SLy4. The EC rates correspond to the stellar 
conditions (a) pY e = 10 7 mol/cm 3 , T g = 3 GK, and (b) pY e = 
10 9 mol/cm 3 , T 9 = 10 GK. 



apply to 62 Ni and 64 Ni. Finally, the rates in Fig. 1231 for 
64 Zn with a very small Qec energy, are sensitive to high 
excitation energies that can contribute especially at high 
p values. 

To quantify the measure of the quality of the various 
calculations for the EC rates, the ratios of the EC rates 
are compared in Fig. [25] for the pf-sheil nuclei studied 
in this work and for two stellar conditions correspond- 
ing to pY e = 10 7 mol/cm 3 , T 9 = 3 GK in Fig. HSJa) 
and to pY e = 10 9 mol/cm 3 , T 9 = 10 GK in Fig. HSJb). 
Because we do not have experimental GT strength distri- 
butions from high resolution charge-exchange reactions, 
(d, 2 He) or (t, 3 He), for all of these nuclei, we show in this 
comparison the relative EC rates using the QRPA-SLy4 
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results as a reference, k>g 10 (A E c/AEC,QR.PA-SL y 4)- As it 
was noticed in Ref. [32|], the results in Fig. 1231 confirm 
that the SM calculations reproduce fairly well the rates 
from the high resolution data, while our QRPA results 
improve significantly the results from QRPA-Mollcr, ex- 
cept in the case of 48 Ti. The comparison of the calcu- 
lations with the rates obtained from (n,p) data shows a 
tendency to underestimate them in most cases. It is also 
worth mentioning the better agreement with the data 
with the stellar conditions given in Fig. 125T b). where 
the EC rates are less sensitive to the details of the GT 
strength distributions. 

The quality of the calculations can be further quanti- 
fied with the help of the average ratios of all nuclei defined 
by 



Zn — > Cu 
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Aec = nY, 1o %w [Ai(th)/Ai(exp)] , 



(13) 



or with the average of the absolute ratios that avoid pos- 
itive and negative cancellations and gives us the average 
order of magnitude of these ratios, 
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These average results are given in Tables III and IV for 
the two stellar conditions discussed above. In these tables 
one can see the ratios of various theoretical EC rates with 
respect to the experimental rates for p/-shell nuclei. We 
consider two sets of experimental rates, in the first row 
of each model the ratios are calculated with respect to 
the rates obtained from the GT strength distributions 
extracted from (n,p) charge- exchange reactions. In the 
second row the ratios are calculated with respect to the 
high resolution charge exchange reactions (h.r.). The last 
two columns correspond to the average of the ratios and 
to the average of the absolute ratios, as defined in Eqs. 
(IT31 and fll~4")) . respectively. In the case of (n,p) data 
the average involves 11 nuclei, whereas in the case (h.r.) 
there are 5 nuclei involved. 

Table III shows that for these stellar conditions, the 
different theoretical models produce comparable |Aec| 
average ratios in the case of (n.p) data, with a little bit 
better agreement for the QRPA calculations with Skyrme 
forces. On the other hand, the best agreement in the case 
of high resolution data is obtained with the SM calcula- 
tions, although the agreement in the QRPA-Skyrme cases 
is also quite good. In the stellar conditions of Table IV 
the quality of the calculations improve and the agreement 
with experiment is much better. The average of the ab- 
solute ratios of the different models with respect to the 
{n, p) data are comparable. On the other hand, the agree- 
ment with the high resolution data is especially good for 
the SM with the GXPF1 force. The agreement obtained 
from SM with KB3G and from QRPA with SLy4 and SG2 
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FIG. 26: (Color online) Electron-capture rates in 64 Zn calcu- 
lated from QRPA-SLy4 at various densities and temperatures. 
The total rates are decomposed into their contributions from 
decays of the ground state + and the 2 + excited state. 



Skyrme forces is comparable, while the QRPA-Moller re- 
sults are somewhat worse. 

As it has been mentioned, ECs from excited states in 
the parent nuclei have not been considered in this work on 
the basis of the typical excitation energies of the lowest 
2 + excited states that are in most cases larger than 1 
MeV and the T-range considered that is not sufficiently 
high to populate those states significantly. 

Nevertheless, in order to check to what extent these 
contributions are small, we have calculated the EC rates 
originated by transitions from the 2+ excited state in the 
case of 64 Zn, where E 2 + = 0.992 MeV [5l| is one of the 
lowest excitation energy in our set of nuclei. We compare 
in Fig. [55] these EC rates with those from the ground 
state in the case of QRPA-SLy4. In this figure we can 
see the various contributions to the EC rates for several 
densities as a function of T coming from the transitions 
from the ground state 0j s — > 1 + and from the transitions 
from the excited state 2 + — > 1 + , 2 + , 3 + in the parent nu- 
cleus, as well as the sum of these contributions (labeled 
as total rate). The figure also shows for comparison the 
EC rates calculated from the ground state alone with a 
probability of population set to one (g.s. only), which 
correspond to the QRPA-SLy4 results in Fig. [23 The 
results show that EC from the ground state is the domi- 
nant contribution except at high T, where contributions 
from excited states are comparable. We also observe that 
the EC from only the ground state is always larger than 
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the contribution 0+ — >• 1 + to the total rate because of 
the depopulation of the ground state when the 2 + state 
becomes populated. 



TABLE II: Electron chemical potentials fi e (MeV) for se- 
lected values of densities pY e (niol/cm 3 ) and temperatures T% 
(GK). 
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0.672 
2.437 
11.116 


0.299 
2.355 
11.098 


0.091 
2.192 
11.063 


0.042 
1.952 
11.011 


0.020 
1.493 
10.898 



enough, the EC rates are sensitive to all the spectrum. 
Then, the whole description of the GT strength distribu- 
tion is more important than a detailed description of the 
low-lying spectrum. Since QRPA reproduces reasonably 
well the global behavior of the GT strength distributions, 
it is expected to be a good approach, especially for high 
p and T conditions. 

(3 + decays and EC from thermally populated excited 
states in the parent nuclei are expected to occur at high 
T. Contributions to the EC rates from the low-lying 
excited state in 64 Zn have been evaluated and compared 
with the contributions from the ground state in the range 
of p and T considered in this work. A systematic study of 
these contributions will be very important for more ex- 
treme stellar conditions, but in this work they can be 
safely neglected. Moreover, the main purpose of this 
work is to compare QRPA results with benchmark SM 
calculations and to EC rates extracted from the measured 
GT strength distributions in the laboratory that do not 
contain those type of contributions. 



IV. CONCLUSIONS AND FINAL REMARKS 

In this work we have evaluated continuum electron- 
capture rates at different density and temperature con- 
ditions holding in stellar scenarios. This study is per- 
formed on a set of p/-shell nuclei representative of the 
constituents in presupernova formations (i.e., Sc, Ti, 
51 V, 54 Fe, 55 Mn, 56 Fe, 58 Ni, 59 Co, 60 Ni, 62 Ni, 64 Ni, 
and 64 Zn). The nuclear structure involved in the cal- 
culation of the energy distribution of the Gamow- Teller 
strength is described within a selfconsistent deformed 
HF+BCS+QRPA formalism with density-dependent ef- 
fective Skyrme interactions and spin-isospin residual in- 
teractions. 

We find that the present QRPA calculation is able to 
reproduce the main features of the GT distributions ex- 
tracted in these nuclei from charge- exchange reactions. 
Comparison of our results for the EC rates with SM cal- 
culations and other QRPA results analyzed in Ref. [32J 
shows that, in general, the agreement with experiment 
is best in the case of SM calculations, but the present 
QRPA calculations with Skyrme forces are of similar 
quality in most cases, clearly improving the agreement 
with experiment with respect to QRPA-Moller. Thus, 
the results in this work provide additional information 
on the performance of QRPA-based models and refine the 
conclusions in Ref. [32|], where it was found that QRPA 
calculations based on the model developed in Ref. [35[ 
produce systematically much larger deviations from the 
data than SM calculations. This work demonstrates that 
this is not the case for all QRPA type calculations. 

We have studied the sensitivity of the EC rates to both 
p and T. At low p (low Fermi energies) and low T (sharp 
shape of the energy distribution of the electrons), the 
rates are very sensitive to details of the GT strength of 
the low-lying excitations and therefore to model calcula- 
tions. On the other hand, when the p and T are high 
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TABLE III: Ratios of various theoretical EC rates with respect to the experimental rates, log 10 (AEC,th/AEC,ex P ), for p/-shell 
nuclei. The experimental EC rates of reference correspond to the GT strength distributions extracted from (n,p) or high 
resolution charge exchange reactions (h.r.). The stellar density and temperature conditions are pY e — 10 7 mol/cm 3 , Tg — 3 
GK. The last two columns correspond to the average of the ratios Aec and to the average of the absolute ratios |Aec|, as 
defined in Eqs. I]13[l and (|14|l . respectively. 



'Sc 



s Ti 



L V 



*Fe 



3 Mn 



'Fe 



s Ki 



J Co 



'Ni 



-Ni 



»Ni 



*Zn A E 



|A E 



SM-KB3G 


'h.r.) 


-4.518 


-2.544 
-0.130 


-1.320 
-0.040 


-1.439 


-2.228 


-2.376 


-1.526 
-0.460 


-1.671 


-0.980 


-1.470 


-0.223 
0.159 


-0.589 


-1.845 1.845 
-0.212 0.276 


SM-GXPF1 


[n,p) 
'h.r.) 


-6.357 


-2.596 
-0.182 


-1.217 
0.063 


-1.490 


-2.081 


-1.869 


-1.746 
-0.680 


-1.924 


-0.298 


-1.082 


-0.334 

0.048 


-0.446 


-1.909 1.909 
-0.239 0.284 


QRPA-Moller 


[n,p) 

'h.r.) 


1.219 


-0.240 

2.174 


0.642 
1.922 


-3.518 


-1.864 


-4.551 


-3.151 
-2.084 


-0.820 


-1.849 


-1.417 


-1.316 
-0.934 


-1.034 


-1.533 1.872 
0.009 1.630 


QRPA-SLy4 


'h.r.) 


0.175 


0.217 
2.631 


-1.066 
0.214 


-0.296 


-3.155 


-3.063 


-0.936 
0.131 


-3.169 


-0.620 


-0.979 


-0.685 
-0.303 


-0.543 


-1.234 1.306 
0.426 0.764 


QRPA-SG2 


Xp) 
'h.r.) 


0.172 


0.211 

2.625 


-1.001 
0.279 


-0.296 


-3.103 


-2.924 


-0.988 
0.078 


-2.826 


-0.010 


-0.799 


-0.590 
-0.208 


-0.041 


-1.105 1.175 
0.547 0.646 



TABLE IV: Same as in Table III, but for the stellar conditions pY e = 10 9 mol/cm 3 , T 9 = 10 GK. 







45 Sc 


48 Ti 


Sly 


54 Fe 


55 Mn 


56 Fe 


58 Ni 


59 Co 


60 Ni 


62 Ni 


64 Ni 


64 Zn 


Aec 






|Aec| 


SM-KB3G 


'n,p 


-0.917 


-0.846 


-0.075 


-0.190 


-0.286 


-0.359 


0.186 


0.061 


-0.197 


-0.142 


0.158 




-0.237 


0.311 




'h.r/ 




-0.138 


0.126 








0.245 








0.147 


-0.597 


-0.043 


0.251 


SM-GXPF1 


n,p 


-0.976 


-0.813 


-0.139 


-0.337 


-0.385 


-0.346 


0.014 


-0.172 


-0.073 


-0.093 


0.047 




-0.298 


0.309 




'h.r/ 




-0.105 


0.062 








0.072 








0.036 


-0.211 


-0.029 


0.097 


QRPA-Moller 


'n,p 


1.050 


-0.125 


0.471 


-0.399 


0.134 


-0.628 


0.023 


0.211 


-0.169 


-0.242 


-0.472 




-0.013 


0.357 




'h.r/ 




0.583 


0.672 








0.082 








-0.483 


-0.039 


0.163 


0.372 


QRPA-SLy4 


n,p 


0.369 


-0.161 


-0.644 


-0.466 


-0.604 


-0.662 


0.040 


-0.401 


0.018 


0.140 


-0.108 




-0.226 


0.328 




'h.r/ 




0.547 


-0.443 








0.099 








-0.118 


-0.206 


-0.024 


0.283 


QRPA-SG2 


n,p 


0.374 


-0.103 


-0.430 


-0.318 


-0.580 


-0.476 


0.138 


-0.253 


0.038 


-0.020 


-0.045 




-0.152 


0.252 




;h.r.; 




0.605 


-0.228 








0.197 








-0.055 


0.073 


0.118 


0.232 



[1] E. M. Burbidge, G. R. Burbidge, W. A. Fowler, and F. [10 

Hoyle, Rev. Mod. Phys. 29, 547 (1957); G. Wallerstein [11 

et al, Rev. Mod. Phys. 69, 995 (1997). 

[2] G. M. Fuller, W. A. Fowler, and M. J. Newman, Ap. J. [12 

Suppl. 42, 447 (1980); Ap. J. 252, 715 (1982); Ap. J. [13 

Suppl. 48, 279 (1982); Ap. J. 293, 1 (1985). 

[3] M. B. Aufderheide, I. Fushiki, S. E. Woosley, and D. H. [14 

Hartmann, Ap. J. Suppl. 91, 389 (1994); A. Heger, S. E. [15 

Woosley, G. Martmez-Pinedo, and K. Langanke, Ap. J. [16 

560, 307 (2001). ' [17 

[4] K. Langanke and G. Martmez-Pinedo, Rev. Mod. Phys. [18 

75, 819 (2003). [19 

[5] W. P. Alford et al, Nucl. Phys. A531, 97 (1991). [20 

[6] K. Yako et al, Phys. Rev. Lett. 103, 012503 (2009). 

[7] S. Rakers et al, Phys. Rev. C 70, 054302 (2004). [21 

[8] W.P. Alford et al, Phys. Rev. C 48, 2818 (1993). [22 

[9] C. Baumer et al, Phys. Rev. C 68, 031303(R) (2003). 



T. Ronnqvist et al, Nucl. Phys. A563, 225 (1993). 

M. C. Vetterli et al, Phys. Rev. Lett. 59, 439 (1987); 

Phys. Rev. C 40 559 (1989). 

S. El-Kateb et al, Phys. Rev. C 49, 3128 (1994). 

M. Hagemann et al, Phys. Rev. C 71, 014606 (2005); 

Phys. Lett. B 579, 251 (2004). 

A. L. Cole et al, Phys. Rev. C 74, 034333 (2006). 

A. L. Williams et al, Phys. Rev. C 51, 1144 (1995). 

N. Anantaraman et al, Phys. Rev. C 78, 065803 (2008). 

L. Popescu et al, Phys. Rev. C 75, 054312 (2007). 

E.-W. Grewe et al, Phys. Rev. C 77, 064303 (2008). 

G. W. Hitt et al, Phys. Rev. C 80, 014313 (2009). 

Y. Fujita, B. Rubio, and W. Gelletly, Prog. Part. Nucl. 

Phys. 66, 549 (2011). 

F. Osterfeld, Rev. Mod. Phys. 64, 491 (1992). 

S. E. Koonin, D. J. Dean, and K. Langanke, Phys. Rep. 

278, 1 (1997). 



18 



[23] K. Langanke and G. Martmez-Pinedo, Nucl. Phys. 

A673, 481 (2000). [41 

[24] K. Langanke and G. Martmez-Pinedo, At. Data Nucl. 

Data Tables 79, 1 (2001). [42 

[25] T. Suzuki et al, Phys. Rev. C 79, 061603(R) (2009); 

Phys. Rev. C 83, 044619 (2011). [43 

[26] P. Moller, J. R. Nix, and K. -L. Kratz, At. Data Nucl. 

Data Tables 66, 131 (1997). [44 

[27] J.-U. Nabi and H. V. Klapdor-Kleingrothaus, At. Data [45 

Nucl. Data Tables 71, 149 (1999); 88, 237 (2004). 
[28] J.-U. Nabi, Eur. J. A. 40, 223 (2009); Astrophys. Space [46 

Sci. 331, 537 (2010); Eur. J. A. 48, 84 (2012). [47 

[29] N. Paar, T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. Mg 

C 69, 054303 (2004); N. Paar, G. Colo, E. Khan, and D. 

Vretenar, Phys. Rev. C 80, 055801 (2009). [49 

[30] P. Sarriguren, Phys. Lett. B 680, 438 (2009); Phys. Rev. r 50 

C 83, 025801 (2011). 
[31] A. F. Fantina, E. Khan, G. Colo, N. Paar, and D. Vrete- [51 

nar, Phys. Rev. C 86, 035805 (2012). 
[32] A. L. Cole et al, Phys. Rev. C 86, 015809 (2012). [ 52 

[33] A. Poves, J. Sanchez- Solano, E. Caurier, and F. Nowacki, 

Nucl. Phys. A694, 157 (2001). [53 

[34] M. Honma, T. Otsuka, B. A. Brown, and T. Mizusaki, 

Phys. Rev. C 65, 061301(R) (2002); 69, 034335 (2004). [54 

[35] J. Krumlinde and P. Moller, Nucl. Phys. A417, 419 

(1984); P. Moller and J. Randrup, Nucl. Phys. A514, r 55 

1 (1990). 
[36] P. Moller, J. R. Nix, W. D. Myers, and W. J. Swiatecki, r 56 

At. Data Nucl. Data Tables 59, 185 (1995). [57 

[37] J. A. Halbleib and R. A. Sorensen, Nucl. Phys. A98, 542 

(1967). [58 

[38] K. Muto, E. Bender, and H.V. Klapdor, Z. Phys. A 

333, 25 (1989); K. Muto, E. Bender, T. Oda, and H.V. [59 

Klapdor-Kleingrothaus, Z. Phys. A 341, 407 (1992). 
[39] F. Frisk, I. Hamamoto, and X. Z. Zhang, Phys. Rev. C rgg' 

52, 2468 (1995). r 61 

[40] P. Sarriguren, E. Moya de Guerra, A. Escuderos, and A. 



C. Carrizo, Nucl. Phys. A635, 55 (1998). 

P. Sarriguren, E. Moya de Guerra, and A. Escuderos, 

Nucl. Phys. A658, 13 (1999). 

P. Sarriguren, E. Moya de Guerra, and A. Escuderos, 

Nucl. Phys. A691, 631 (2001). 

P. Sarriguren, E. Moya de Guerra, and A. Escuderos, 

Phys. Rev. C 64, 064306 (2001). 

I. N. Borzov, Nucl. Phys. A 777, 645 (2006). 

P. Sarriguren, E. Moya de Guerra, R. Alvarez- Rodriguez, 

Nucl. Phys. A716, 230 (2003). 

E. Poirier et al, Phys. Rev. C 69, 034307 (2004). 

E. Nacher et al, Phys. Rev. Lett. 92, 232501 (2004). 

P. Sarriguren, R. Alvarez-Rodriguez, and E. Moya de 

Guerra, Eur. Phys. J. A 24, 193 (2005). 

P. Sarriguren, Phys. Rev. C 79, 044315 (2009). 

P. Sarriguren and J. Pereira, Phys. Rev. C 81, 064314 

(2010). 

Evaluated Nuclear Structure Data File (ENSDF), 

http://www.nndc.bnl.gov/ensdf/ 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. 

Schaeffer, Nucl. Phys. A635, 231 (1998). 

N. Van Giai and H. Sagawa, Phys. Lett. B 106, 379 

(1981). 

D. Vautherin and D. M. Brink, Phys. Rev. C 5, 626 
(1972); D. Vautherin, Phys. Rev. C 7, 296 (1973). 

H. Flocard, P. Quentin, A. K. Kerman, and D. Vautherin, 

Nucl. Phys. A203, 433 (1973). 

H. Homma et al, Phys. Rev. C 54, 2972 (1996). 

A. Bohr and B. Mottelson, Nuclear Structure, (Benjamin, 

New York 1975). 

P.B. Radha, D.J. Dean, S.E. Koonin, K. Langanke, and 

P. Vogel, Phys. Rev. C 56, 3079 (1997). 

E. Caurier, K. Langanke, G. Martmez-Pinedo, and F. 
Nowacki, Nucl. Phys. A653, 439 (1999). 

W. R. Hix et al, Phys. Rev. Lett. 91, 201102 (2003). 
K. Iwamoto et al, Ap. J. Suppl. 125, 439 (1999). 



